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Abstract. A phenomenological power spectrum of primor- 
dial density perturbations has been constructed by using both 
COBE data to probe the large wavelength region, and a dou- 
ble power law, locally deduced from galaxy catalogs, which 
describes the matter correlation function up to tens of Mega- 
parsec. The shape of the spectrum P{k) of density fluctuations 
exhibits a peak that singles out a characteristic wavelength 
Apeak proportional to the cutoff radius Ro in the matter distri- 
bution (comparable to the distance at which matter becomes 
anticorrelated) . From a least squares fit to COBE's angular 
correlation function we obtain i?o = 35 ± 12/1" ^Mpc for the 
correlation length, and n = 0.76 ± 0.3 for the spectral index 
of P{k) in the large wavelength region. The inferred scale in 
the spectrum is Apeak = 51 ± 18/i~^Mpc. This number agrees 
with that derived from the analysis of the correlation function 
of matter and with a preferred scale identified in IRAS PSC. 

Key words: Galaxies: clustering - Cosmology: cosmic mi- 
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. ^ 1. Introduction 



After the positive detection of anisotropies in the cosmic mi- 
crowave background radiation (CMB) by COBE-DMR (Ben- 
nett et al. 1992; Smoot et al. 1992), it is now possible to directly 
probe the spectrum of primordial density fiuctuations (PDF). 
Early attempts to place restrictions on the spectrum of PDF re- 
lied mainly on the existing data on large scale matter distribu- 
tion or on theoretical prejudices, depending on what region of 
the spectrum was explored. While matter distribution is useful 
in determining the shape of the 'evolved' spectrum at relatively 
large wave-numbers, large angular scale CMB anisotropies are 
sensitive to the primordial spectrum at low wave-numbers. 

Analysis of the QDOT IRAS surveys and the angular two- 
point autocorrelation function w{6) (Peacock 1991) shows that 
both data sets are consistent with a spectral function P{k) that 
approaches A;~^'^ for large k. However, increasing the depth of 
the surveys has revealed the presence of a break in ^(r) (the 
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two-point galaxy-galaxy correlation function; Calzetti et al. 
1992), the existence of a region where ^(r) is negative (Guzzo 
et al. 1991), and also the possible detection of a cutoff radius 
_Ro of order 30 Mpc. In this paper we make use of the 
double power law in 1 -|- ^(r) to derive a properly normalized 
power spectrum for large k, and we match it to a power law 
P{k) oc fc" for small k. The matching point is selected so that 
the power spectrum makes a smooth transition from the pri- 
mordial power law to the 'evolved' spectrum. As a result, the 
peak in the spectrum turns out to be of the same order of 
magnitude as Ro . The cutoff radius in the matter distribution 
i?o and the spectral index n are then the only two important 
quantities, and we can treat them as free parameters. A cal- 
culation of anisotropies in the CMB with this model for the 
power spectrum allows us to fit these parameters to COBE's 
angular correlation function. The quadrupole amplitude allows 
for an independent check for consistency. 

The fitting turns out to be very satisfactory and provides 
for Ro a best value of 35 ± 12/i~^Mpc very close to the optimal 
value inferred from observations of the galaxy-galaxy corre- 
lation function (Calzetti et al. 1991). The best value for the 
primordial spectral index is found to be 0.76 ± 0.3, with a 
central value somewhat lower than predicted by standard ex- 
ponential inflation but, because our error bars, still consistent 
with it and with Smoot et al. (1992). A relatively low index is 
not peculiar to the model because for a pure power-law spec- 
trum we find similar results. With respect to this, we notice 
that an attempt to fit COBE data to a power law spectrum of 
density fluctuations within CDM models results in a spectral 
index 0.6 < n < 1.0 (Adams, et al. 1992). Kashlinsky (1992) 
also finds n < 1.0 in order to satisfy data on both the large 
scale galaxy distribution and large scale streaming motions. 

A feature which, on the other hand, is peculiar to our 
model is the close relationship between the cutoff radius of 
the matter correlation function and the peak in the spec- 
trum. As a matter of fact, the maximum of P{k) is located 
at Apeak = 51± 18/i~^Mpc. This number agrees with the privi- 
leged length scale of about 50/1"^^ Mpc found in the IRAS Point 
Source Catalog by Fabbri and Natale (1993a). A further anal- 
ysis (Fabbri and Natale 1993b) provides evidence for the in- 
terpretation of this scale as a peak in the product k^^'^Pik). 
However, it does not agree with the very large scales recently 
inferred from redshift surveys (Fisher et al. 1993; Einasto et 
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al. 1993; Jing & Valdarnini 1993). This point will be discussed 
in Sect. 5. 



2. Power spectrum 

To construct a phcnomcnological spectrum wo started from the 
autocorrelation of the density fluctuation Sp/p, which is equal 
(up to a delta function at the origin) to the particle correlation 
function ^(r) and can be parametrized as follows (Calzetti et 
al. 1992): 
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To make the calculations independent of Hubble's constant 
Hq = lOO/i Km s~^ Mpc^^ we adopt the change of variable 
X = hr. The correlation function admits four independent pa- 
rameters, which we choose to be the spectral indexes 71 and 72 
and the transition points Xt and xq. The amplitudes Ai and A2 
are then calculated in terms of these parameters by demanding 
continuity of 1 + ^{x) ai x = Xt and the normalization condi- 
tion on the correlation function ^(r)r^dr = 0. The resulting 
formulae for the amplitudes are: 
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Three parameters are determined by the existing data, 
namely 71 = 1.8, 72 = 0.8, xt = 3.0 Mpc. It turns out that 
our results are quite insensitive to the precise values of 71 and 

Xt, and moderate differences are found varying 72 in the range 
0.7 — 0.9. On the other hand, the value of the cutoff radius 
Xq — hRa is important for our calculations because the depen- 
dency of A\ and A2 on this parameter affects the normalization 
of P{k) and thereby the CMB anisotropies. We treat xo as a 
free parameter since it is still poorly determined. 

Let us now define the spectrum of PDF in terms of the rms 
density fluctuations: 



P{k)^l(ri)d?k, 



where A; is a dimensionless wavenumber defined as: 
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with k the physical comoving wavenumber, a = 
1/6000 Mpc~^, and rj the conformal time. The temporal evolu- 
tion of perturbations is described by ^^(t)), which in the linear 
regime is independent of wavenumber, and for fl = 1 can be set 
equal to if /IQ with rjo = 1. The spectrum of PDFs, being the 
Fourier transform of the matter correlation function (Peebles 
1980), results in 
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The above integration was done numerically using a combi- 
nation of both the Clcnshaw-Curtis and Gauss-Kronrod meth- 
ods which are adequate for an oscillating integrand (Piessens et 
al. 1983). The results of the Fourier trausform shows the pres- 
ence of a peak in P{k) at a wavenumber fcpeak inversely propor- 
tional to Xq^ (/Cpeak ~ 2.56 X 10''Mpc/a;o) . Beyond the peak the 
function P{k) shows an oscillatory behavior (due to the sharp 
edge cutoff of £,(r)) at first, but eventually reaches a power 
law behavior with a negative slope P{k) oc k~^''^, in agreement 
with the results of Peacock (1991). For small wavenumbers 
k < fepeak we find P{k) oc k^. This asymptotic law would be 
found from any correlation function which is forced to vanish 
beyond a maximum radius I'o, or decreases as x ' with 7 > 5, 
provided the integral of ^ over space vanishes. This behavior 
does not agree with COBE's data. However, since the domain 
where the galaxy-galaxy correlation function appreciably dif- 
fers from zero is limited to relatively small scales, P{k) can 
be modified in the small-fc region without appreciably affect- 
ing ^(x). Therefore we replaced the P{k) cx k^ law by a power 
function with a spectral index n left as a free parameter to 
fit with COBE's data. In order to minimize the perturbation 
to ^(a;) we extrapolated the spectrum in Eq. (6) down to the 
wavenumber km where its slope assumes the desired value of 
n. For k < km the power law function P{k) = Pm{k/km)" 
was adopted. The parameter Pm is the value of the original, 
5— derived spectrum in A; = km- Continuity of both the spectral 
function and its first derivative is therefore guaranteed. Figure 
1 illustrates the form of the function P{k) used here for our 
best values of 2:0, n and for 72 — 0.8. We have found a power 
series solution to the Fourier transform in Eq. (6) neglecting 
the region where ^(r) has slope —1.8: 
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which can be conveniently used to find the matching point km 
of the function an its first derivative. An approximation which 
is sufficient to our purpose is 
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Since Eqs. (6)-(8) do not contain Xt and 71, they cannot 
describe the large-fe region. However, this region is quite inef- 
fective to infiuence the CMB anisotropies. 

3. CMB Anisotropies 

With a properly normalized spectrum of PDF, one can cal- 
culate the CMB large scale fiuctuations and angular correla- 
tion function. The anisotropy field is conveniently expressed in 
terms of the coefficients in a harmonic expansion: 
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of P{k) goes from fcpeak = 1-3 x 10^ to 256. In our calculations 
the shape of P{k) around the maximum is not so important for 
the determination of the anisotropies of the CMB, other than 
providing an overall normalization. 

4. Results 

The possibility of using the DMR angular correlation data 
to find restrictions on the spectrum of PDF was exploited 
here by fitting the theoretical angular correlation function 
C{6,3.1°;xo,n) without any arbitrary normalization factors to 
that obtained by COBE-DMR. For each value of 72 the pa- 
rameters resulting from the fit are the cutofi' radius xo and the 
spectral index n of the long-wave side of PDF spectrum. 

A numerical minimization of the function using CERN's 
MINUIT package (James & Roos 1975) gives the results re- 
ported in Table 1. 



Fig. 1. Spectrum of primordial density fluctuations for xo 
Mpc, n = 0.76 and 72 = 0.8 
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Table 1. Best fitted parameters 
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with the dipole terms = 1) excluded. For random per- 
turbation fields the variance of aim defines the angular power 
spectrum: 



(2£+l)(oL>- 
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Experiments such as COBE-DMR which produce a sky 
map of the cosmic background radiation give a direct mea- 
sure of the angular correlation function C{6,a), which can be 
theoretically predicted by 



(13) 



where Pt(cos9) are the Legendre polynomials, and the 'fil- 
ter' function for a Gaussian beamwidth cr = l/(£o + |) radians 
is 



Wg. = exp 



2(4 + 



(14) 



For COBE-DMR filtering efi'ectively starts at io « 18 (cor- 
responding to an angular resolution of 3.1°), but in our case 
terms up to £ fs 60 should be included for an accurate calcular 
tion. 

The imprint left on the CMB produced by metric pertur- 
bations during the epoch of last scattering, (Sachs & Wolfe, 
1967) can be calculated by evaluating the integral: 



«?=(f)'(2£+l) j P{k) 



k{k) 



dk, 
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with je the spherical Bessel functions. For a cutoff radius xo 
lying in the range from 20 to 100 Mpc the location of the peak 
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36.4 ± 12 


0.76 ± 0.3 


64.055 


0.8 


34.8 ± 12 


0.76 ± 0.3 


64.055 


0.9 


33.5 ± 12 


0.76 ± 0.3 


64.055 



Fig. 2. Angular power spectrum for the best fit values as in Fig. 1 

The fitted parameters are found to be fairly insensitive to 
the precise value of 72, so that we can safely assume that xo = 
35±12 Mpc and n = 0.76±0.3. The errors quoted correspond to 
the 68% confidence level. The reported are for 69 — 2 degrees 
of freedom. Figures 1, 2 and 3 respectively show the power of 
density fiuctuations, the CMB angular power spectrum, and 
the correlation function obtained by the fit for the 72 = 0.8 
case. In Fig. 3 the DMR data is also included in the plot. 
Figures 4 and 5 give contour plots of near the minimum. 
The contours show that the two-parameter fit exhibits a strong 
correlation. For 72 = 0.8 we have xo = 0.48-1-45. 54n Mpc, with 
r = 0.9996. This effect is built-in by the way we normalize the 
galaxy-galaxy correlation function. From Eqs. (2) and (3) it 
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Fig. 3. Angular correlation function for best fit values as in Fig. 
1 and COBE-DMR data points. The dipole and quadrupole terms 
have been removed 



Fig. 5. Contour plots for a) 72 = 0.7, and b) 0.9. The reported 
contours are for Ax^ = 2.3 and 11.8 



appears that increasing xq the amplitude of the galaxy-galaxy 
correlation function and thereby of P{k) increases, so that we 
need to increase n in order to suppress the excess of long-wave 
power. 



a free parameter of the fit (which in the end is used to ronor- 
malize the spectrum). The results arc — 64.06, quadrupole 
amplitude, Qrms-ps = (AT)2 = ro(ai/47r)^/^ of 16±8 /xK and 
n = 0.8 ± 0.3. If one fixes n to 1.0 then increases to 64.88 
while the quadrupole amplitude decreases to Q = 15 ± 5 fjK. 

We have checked the consistency of our results with the 
independent measurement of the RMS quadrupole anisotropy 
by COBEl-DMR. As can be seen in Fig. 6, the quadrupole 
evaluated at the optimal values of and n agrees with COBE's 

Qrras — ps • 



Fig. 4. Contour plot of x^ around the minimum for 72 = 0.8. The 
levels are for Ax^ = 2.3,6.17 and 11.8, corresponding to 68.3%, 
95.4% and 99.7% 

One could argue that if the cutoff radius xo is an indepen- 
dently well known parameter, then it should be fixed. Taking 
Xq = 30 Mpc (the best value of Calzetti et al. 1992) the spec- 
tral index obtained from the fit is more tightly constrained: 
n = 0.64 ± 0.02 with a of 64.2. This is far from unity. 

Alternatively one could fix n to 1.0 and ask for the best xo 
that fits the data, still satisfying the second of the conditions 
above. This results in xo = 44 ± l/i~^ Mpc with = 64.88. 
This radius implies Apeak = 65/i~^ Mpc. 

To compare with the numbers obtained by the COBE- 
DMR group we have made the fit using a 'pure', unnormalized 
power law P{k) oc k", and left the quadrupole amplitude ol as 



Fig. 6. Quadrupole amplitudes derived from the model with 
72 = 0.8 for various values of n. Solid lines are curves of fixed n, 
with n = 0.6, 0.7, 0.8, 0.9, 1.0, 1.2 and 1.-1 going from the steepest 
to the least steep curve. The quadrupole amplitudes corresponding 
to the optimal value of xq axe indicated by squares. The error bars 
at the absolute minimum correspond to the 68% confidence level. 
The horizontal lines indicate the 68% confidence region of COBE's 

Qrms — ps 
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5. Discussion 

The cutoff radius Jio = 35 ± Uh'^ Mpc ag rees very well with 
the result found by Calzetti et al. (1992) on the Zwicky/CfAl, 
SSRS, and IRAS PSC catalogs. Since _Ro corresponds to a 
wavelength of 51h^ Mpc it agrees with the preferred wave- 
length found by Fabbri and Natale (1993a) in IRAS PSC, 
A = 50h~^ Mpc. However, the best value for the primordial 
spectral index n is found to be somewhat lower than the best 
value declared by Smoot et al. (1992), and from that expected 
from the scenario of standard inflation. The possibility of lower 
values of n from DMR data has been noticed also by Adams 
et al. (1992) and Kashlinsky (1992). 

Is our analysis consistent with standard inflation? Inspec- 
tion of our contours in Figs. 4 and 5 shows that it is so only 
paying the price of a rather large value of the cutoff radius, 
Ro w 45/i~^ Mpc. Such a value weakens the agreement with 
Calzetti et al. (1992). However we should notice that our model 
- obtained by the requirement of a minimum perturbation of 
^(r) - forces the test to be more stringent than other tests 
on DMR data since there is no normalization factor as a free 
parameter. We cannot exclude that different spectral shapes 
might provide results more favourable towards standard infla- 
tion, also a pure power-law model with a best-fitted normal- 
ization factor does favor values of n smaller than unity. It is 
worth noticing that so small values of n also arise theoretically 
from power-law inflation (Fabbri et al. 1987). 

While the best value found for the correlation scale is com- 
pletely satisfactory, the peak wavelength found for P{k) is not 
consistent with spectra recently derived from redsliift surveys 
on IRAS galaxies (Fisher et al. 1993), clusters (Einasto et al. 
1993) and both (Jing & Valdarnini 1993). There the power 
spectra appear to increase towards long wavelenths, and a 
turnover appears only at the boundaries of the samples A ~ 150 
or 180/i^^ Mpc. In our view, the significance of such results is 
weakened by the criticism of Bahcall et al. (1993), who show 
that very large distortions must be generated by large-scale 
peculiar velocities. The ratio of the values of the power spec- 
trum at A > 100/i^^ and A < Whr^ Mpc, when calculated in 
redshift space, may be ^ 10 times higher than the 'true' ra- 
tio in real space ! Unfortunately the recovery of the genuine, 
real-space P{k) is model-dependent there. The role of peculiar 
velocities is not significant in our work, as well as in Fabbri 
and Natale (1993a, b). Further, it is not clear whether spec- 
tra peaked around 150h-^ Mpc can be made consistent with 
COBE's data as well as other anisotropy data available around 
1°: Kashlinsky's (1992) results seem to imply that compatibil- 
ity cannot be achieved for so broad spectra, while the goal 
could be reached for sharp but unphysical peaks. The present 
paper shows that the peak wavelength must be much smaller if 
we wish (i) to make a really satisfactory fitting of COBE's data 
and (ii) to minimally perturb the spectrum directly originated 
by ^{r) in the region where it can be measured. 
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